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Abstract. We propose a phenomenological model in which a non-minimal coupling between gravity 

i i and dark matter is present in order to address some of the apparent small scales issues of ACDM model. 

When described in a frame in which gravity dynamics is given by the standard Einstein-Hilbert action, 

the non-minimal coupling translates into an effective pressure for the dark matter component. We 

consider some phenomenological examples and describe both background and linear perturbations. 

(ys We show that the presence of an effective pressure may lead these scenarios to differ from ACDM at 

£■ — , the scales where the non-minimal coupling (and therefore the pressure) is active. In particular two 

l/~) effects are present: a pressure term for the dark matter component that is able to reduce the growth of 

£fs structures at galactic scales, possibly reconciling simulations and observations; an effective interaction 

term between dark matter and baryons that could explain observed correlations between the two 

OA components of the cosmic fluid within Tully-Fisher analysis. 
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1 Introduction 

Our present understanding of the universe is encoded in the so called ACDM model which is able to 
explain in a simple and elegant way the evolution and the formation of structures. However, despite 
its great successes [1] and successful tests at the solar system scale [2], General Relativity is still 
far from being the ultimate theory of our universe. In particular, the present day energy density is 
dominated by dark fluid(s) whose nature has not been completely understood so far. At the largest 
scales the accelerated expansion is still a puzzle to be solved while, at smaller scales, galaxy clusters 
and galaxies are dominated by a dark matter component whose particle nature is still unknown. 

Dark matter, introduced in the realm of astrophysics in 1934 by Fritz Zwicky [3], is still missing 
both a definitive detection and a unique theoretical candidate. Over the years many particles have 
been proposed as viable candidates for dark matter (see [4-6] for comprehensive presentations) but 
up to now no convincing indirect [7-9] nor direct [10-14] detection has led to a resolutive answer. 
At cosmological and astrophysical scales there is no need for a particle description of dark matter 
and its dynamics is described via a fluid approximation. This fluid is obtained from the particle 
description through some averaging procedure that keeps track only of those microscopic properties 
that are relevant from a macroscopic point of view, like pressure, velocity and density. 

The most accredited model, the Cold Dark Matter paradigm (CDM), describes dark matter as a 
non-interacting, pressureless fluid with small dispersion velocity. As already pointed out, despite its 
simplicity this model is able to describe the formation and evolution of structures with remarkably 
high precision. However, although its success, recent observations seem to indicate that the CDM 
picture may not be able, as it stands, to explain in a natural way dark matter dynamics at clusters and 
galactic scales. In particular some tension exists between dark matter simulations and observations, 
with regard to both the density profiles of dark matter halos [15-19] and for the number of predicted 
substructures inside a given host halo [20-23]. 1 Indeed, many attempts have been made to solve these 
issues: baryons feedback on dark matter distribution [25-28], self interacting dark matter [29-31], 
fuzzy cold dark matter [34], frustrated cold dark matter [35] in the context of unified dark matter 
dark energy models or warm dark matter (WDM) [32, 33, 36-39], whose free streaming length is 
able to erase structure formation below it. These problems may be related to the lack of a minimum 
length scale for the clustering of dark matter that hence keeps collapsing down to small scales. In 
other words, the net effect can be interpreted as the insurgence of an effective pressure term in the 
dark matter fluid. Recently it has been proposed a method [40] to test the DM equation of state at 



lr The possibility that these issues may be related to a too low resolution of the simulations has been investigated too 

[24]. 



small scales in a rather model independent way by comparing the gravitational potentials measured 
from galaxy rotational curves and from weak lensing which has been applied to data in [41]. The 
results show a value for the equation of state parameter compatible with zero but with a preferred 
value of w = —1/3 at lcr. From cosmological observations, and hence large scales, the equation of 
state for DM is compatible with zero [42, 43]. Interestingly, even if no incompatibility with w=0 has 
been found, a model with w ^ is still compatible with CMB data [43] . 

Another source of tension concerns empirical tests like the baryonic Tully-Fisher relation [44, 45] 
or the constant galactic surface density [46]: they both hint towards a connection between dark 
matter and baryonic properties; these observations are hardly explained within the ACDM paradigm 
unless the amounts of dark matter and baryons are accurately balanced. On the other hand, dark 
matter-baryon interactions are strongly constrained in the context of particle physics so that it is 
very unlikely they can address these issues. 

Curiosly, a model capable of predicting at least some of these empirical tests is the well known 
MOND paradigm [47, 48] where the mass discrepancy problem is explained via a modified Newtonian 
dynamics without requiring dark matter. MOND is a non relativistic model and several generalizations 
to a fully relativistic theory have been attempted [49-52]. At present, none of them seems however 
able to reproduce with accuracy and most of all, simultaneously, all well established cosmological 
features, like the CMB acoustic peaks, the matter power spectrum, the Bullet cluster (see for example 
[53, 54] for attempts in this direction). 

An alternative appealing solution to CDM problems that may reconcile ACDM successes at large 
scales with MONDian behavior at galactic scales, comes from a new class of non-standard dark matter 
gravity couplings [55-57]. In this scenario the fluid content of the ACDM paradigm is maintained but it 
is rather the dynamics of the dark matter fluid at suitable late times and small scales to be changed. In 
particular, we shall here entail the idea that dark matter could, at suitable scales, have undergone some 
sort of phase transition (e.g. something analogue to a Bosc-Einstcin condensation) and consequently 
developed a finite coherence length of size comparable to the local spacetime curvature. In this sense 
dark matter would be capable to probe gravity - e.g. to "feel" second derivatives of the metric - hence 
becoming non-minimally coupled to curvature terms. 

In this paper we further investigate the model proposed in [57] deriving for the first time its 
cosmological implications. We show how this class of models gives rise to two effects: a pressure term 
for dark matter, able to reduce the growth of structures at small scales, plus an effective interaction 
term between dark matter and baryons that might be able to explain observed correlations between 
the two components. 

The paper is organized as follows. In section 2 we introduce the model, in section 3 we describe 
the field/fluid formalism used and write the equations in their general form. In section 4 we specialize 
them to a flat FRLW universe and study the modifications of our model with respect to ACDM while 
in section 5 numerical results are presented. Finally in section 6 we draw out our conclusions. 

2 Non— minimally coupled dark matter 

The most general-scalar tensor theory that gives second order field equations is encoded in the Horn- 
deski action [58] . Unfortunately due to its generality, it is very hard to build a constrained and reliable 
cosmology using this action (see however [59, 60] for attempts in this sense). In [57] a minimal recipe 
was used to build an action in which the dark matter fluid is non-minimally coupled to curvature 
terms and the astrophysical consequences of this couplings were addressed. The action reads: 

c 3 
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where g is the metric in the Jordan frame (we will use the tilde to indicate quantities in the Jordan 
frame), c is the velocity of light, Gjq is the Newton constant, R and R^ u are the curvature scalar and 



Ricci tensor respectively in the metric g, ascai and aRi c are constants, Pdm is the Dark Matter (DM) 
energy density, ps is the baryonic energy density. T(pdm) and 'E.(pdm) are two generic functions of 
the DM density. Here ascai introduces a coupling between pdm and the curvature; in addition, a^n c 
couples the metric and its derivatives also to the DM four velocity Udm- This model gives rather 
precise predictions in the Newtonian limit: the Newtonian gravitational potential \I/ is sourced not 
only by the energy density but also by its gradient and Laplacian: 

V 2 * = 4ttG n (pdm - ?V 2 S( PM/ ) + a Sca iV 2 T(p DM )) . (2.2) 



These features are interesting ingredients in the tentative solution of ACDM small scales issues as 
they generate corrections that may reduce the cuspiness of the DM density profiles and reduce the 
number of satellite halos. We note that our model was developed in a context which is independent 
from Horndeski; it can still be seen as a subclass of those theories, when only fluid quantities related 
to density and quadrivelocity and their coupling to curvature terms are taken into account. 

We will now further extend the analysis of this class of models, investigating their cosmological 
consequences. In order to do this we choose to move to a frame in which the action for gravity 
is described by the standard Hilbert-Einstein action. When in the Jordan frame the non-minimal 
coupling is taken to be with the Einstein tensor, ascai = l/2aRi c = e, this can be achieved with the 
following metric transformation: 

9tiu = 9nv + V" hp, v = (A(p DM ) 2 - 1) g^ v + B{p DM )u ll u v , (2.3) 

where A(pdm) and B(pdm) are generic functions of the DM density; this is the most general relation 
between metrics that respect both causality and the equivalence principle (WEP) and is called disfor- 
mal transformation [61]. We stress that here h^ v is not only a function of the metric but depends also 
on the DM density itself. If we re-express action (2.1) in terms of the metric g^ v and then expand in 
powers of h^ v up to order 0(h 2 ) we obtain the following action: 

S = She\$\ + Sdm [g, Pdm] + S B [g, Pb] + S int [g, Pb,Pdm], (2.4) 

where She\s] is the standard Einstein-Hilbert gravitational action, SoAilg, Pdm] and Szjm [<7? Adm] 
are the DM and baryonic actions in the metric g, Si n t is 



Sint = ■ - h / ^x^J-g (T£' M + Tg u ) h^ig, pdm), (2.5) 
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and represents a new interaction term which in general involves the metric, dark matter and baryons. 
As usual the stress energy tensor 7j,„ for the i component (DM or baryons) is given by 

v-g "g^v 

In this new frame the effects of the non-minimal coupling have been transferred into a coupling term 
for the stress energy tensors of both dark matter and baryons and h^ which, as stressed before, 
is itself a function of the metric as well as of dark matter fluid variables. This translates directly 
into a coupling between DM and baryons and a self coupling for dark matter. Notice that those 
two couplings are not to be intended in the particle physics sense but they rather emerge from a 
geometrical coupling between dark matter and gravity. 

We further stress here that the action (2.4) is obtained from (2.1) through an expansion and 
hence the two actions are equivalent up to order O(h). This doesn't matter as we could have started 
directly from the action (2.4). 

As argued in the introduction we assume the coupling to be active at late times and small scales 
when dark matter dynamics changes giving rise to a coherence length that activates the non-minimal 
coupling. 



3 Dark matter with an effective pressure 

As a first investigation of this model we neglect the baryon contribution and set the value of the 
conformal function in eq. (2.3) to be A = 1. We also switch to field formalism using standard 
conventions [62]: 

p^p(X,<p), p->p(X,<p), u IM = V IM <p/V2X (3.1) 

where now tp is the dark matter scalar field and X = —g^V ^ipV v (p/2 is its kinetic part. The field 
variables are not fundamental fields, but rather stand for a different representation of the fluid. In 
this sense we must ensure that in absence of a coupling we recover the equations for the standard 
dark matter fluid. In general a scalar field will not behave as a pressureless dust; thus, we have to 
enforce this requirement. In order to do this we introduce a Lagrangian multiplier in the Lagrangian 
for the dark matter field 

C = \{x») (x-~V 2 (<p) ) (3.2) 



where V(ip) is a generic potential so that when we take the variation of the action with respect to A 
we get 

^ = 0^X-iv-M = (3.3) 

which exactly sets the pressure to zero [63, 64]. We conjecture that in order for the interaction to be 
active, the temperature of the cosmological bath (i.e. time or redshift) must be below a certain thresh- 
old. This explains the time dependence of the coupling. However, even after the critical temperature 
has been reached, the density may be too low for the interaction to be efficient. Hence one needs 
high densities in order to make the interaction relevant. In order to implement phenomenologically 
the wished scale and time dependences of the coupling, we fix the parametrization of the coupling 
function in the following way: 

h^(X, ip) = F{ - X, ^ V^V v ip. (3.4) 

The function F gives the scale dependence of the interaction by suppressing it until the density is 
large enough to overcome a given threshold fixed by p*. 
The action (2.4) can then be rewritten as: 

S = S HE [g] + S DM [g, ip, A] + ej d 4 x^£ nmc (X, ip), (3.5) 

where now Sdm depends on {ip, A} rather than on pdm and 

C nmc = _ g^) X ( X + V ( A (3 6) 

P* V 2 / 

is the interaction term generated by the non-minimal coupling. The coupling e is a switch that 
is chosen to be zero at t < t c , t c being the time when the coupling is activated; for t > t c , e 
becomes different from zero, reaching a constant value of about e < 1, as required by the expansion 
of huv to order (e 2 ). In other words, e becoming different from zero indicates the onset of the non- 
minimally coupled epoch. At the same time, for the purpose of the present analysis, which is limited 
to order C(e 2 ), its value must be small compared to one. 2 From a cosmological point of view the 
time dependence of the coupling is required because we want to study how the non-minimal coupling 
modifies the ACDM behavior at small scales and late times only, as we stressed in the introduction. 
To be fully rigorous we should have given e a spatial dependence as well. That would cause the 
activation of the non-minimal coupling to happen at different times and in different regions. Such 
a spatial dependence would give rise to a most interesting phenomenology but strongly dependent 



2 The parameter t is just a phenomenological parameter that should be given dynamically by the non— minimal 
coupling mechanism. In this sense we are not introducing a non-dynamical field and background independence is 
preserved. 



on which powering mechanism is chosen for such a dependence on space. We therefore feel that 
this treatment would go beyond the scope of the present paper, requiring a dedicated analysis in a 
separated work. 

The variation with respect to ip and A gives the two combined equations of motion which together 
specify the fluid dynamics: 

A = V-\ V ) [V{p)p, v - (p + e P )6] , (3.7) 

<p = -V(<p), (3.8) 



where 



= ^V M , s V^u", (3.9) 



and where p^ v is the derivative of the density with respect to the field p. The energy density p and 
pressure p are derived from the total dark matter stress energy tensor: 

V = (A + eCJ c ) V 2 {p)u^u v - eC nmc g^ (3.10) 

by a direct confrontation with the form ad the perfect fluid stress energy tensor. We can identify the 
following thermodynamic quantities: 

p = (A + eCy c ) V 2 {ip) - eC nmc , (3.11) 

p = C nmc , (3.12) 

?/' = V- l {<p)V»<p. (3.13) 

Notice that these quantities are not the same appearing in action (2.4) because here the interaction 
term is directly involved in the definition of both density and pressure. In other terms in (2.4) Pdm 
is the DM density for a pressureless fluid which has some non trivial self interaction while here the 
interaction has been absorbed into the definition of the DM energy, pdm in (2-4) is related to the 
{A, tp} variables by the relation pdm = AF(<yj) 2 . Since our knowledge of the DM distribution comes 
through its gravitational effects, definition (3.11) gives the actual measured DM density. 

We also notice that the pressure defined in (3.12) acts only along the direction defined by the 
fluid four velocity, as can be seen from the fact that the four acceleration a M = u"V ' v u^ is identically 
zero. As a consequence we don't have anisotropic stresses that in more general situations may be 
present [57]. 

The effect of the non-minimal coupling is twofold: on the one side it modifies the dark matter 
energy density while on the other side it introduces a pressure term that would be absent in the 
standard ACDM scenario. We expect both these terms to have relevant cosmological consequences at 
the time and scales of interest, as will be shown below. We anticipate that this model is able, with 
appropriate potential shapes, to reduce the source of the gravitational potential with the consequence 
of smoothing out the overdensities at small scales. The equations for a pressureless dust (3.7)-(3.8) 
in the limit of vanishing coupling are independent on the choice for the potential [63] that is to say 
the density scales like a~ 3 no matter what potential is chosen. This does not represent a problem 
since once the mechanism that generates the non-minimal coupling is fixed, so is the interaction 
term which in turns fixes the shape of the potential. We will not investigate the ultimate nature of 
the non-minimal coupling here but we will give some examples of potentials that could lead to an 
interesting phenomenology. 

4 Cosmological framework: background and linear perturbations 

We will now study the cosmological consequences of our model. We will assume a flat FRLW universe 
filled with the dark matter field plus a cosmological constant. 3 In terms of the conformal time 



3 Since we are interested in cluster/galaxy scales we are not concerned here about the nature of dark energy. We 
thus make the minimal choice of a cosmological constant. 
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Pressure 
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V( V ) = y/ftf M e- KV 


p(z) = 


P(z) = -(P 
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1 
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(ep$ M /p*) = 1 
{epE M / Pt ) = -10- 3 


z c = 5 
z c = 5 



Table 1: Functions and parameters. In the first column are reported the two potentials used in the 
paper while in the second column are reported the related pressures as a function of the redshift. The 
third column reports the used value for the scale function F. The value is set to a constant meaning 
that no scale dependence is present. The last two columns report the value for the effective coupling 
constant ep$ M / p* and the redshift at which the non-minimal coupling is activated. 



dr = dt/a we have that: 



n 2 = ^ [(A + e£T c ) V 2 {^) - eC nmc + p A ] 

A' = -a(T)V- 2 (<p) [V(<p)p, v + m(p + ep)} , 
<ff = a(r)V(ip), 



(4.1) 

(4.2) 
(4.3) 



where primes indicate derivatives with respect to the conformal time and p\ = 3A/(8ttG). The 
equation for A can be rewritten in terms of the more physical quantity p defined by equation (3.11) 
which is what appears on the right hand side of the Friedman equation 4.1 with no cosmological 
constant. In this respect we are defining the density of the dark matter fluid as the quantity which 
plays the role of gravitational source. With this consideration the previous system of equations is: 



n 2 



8nGa 2 



[P + Pa] 



p' + 3H(l + ew(<p))p = 

if 1 = a(r)V(ip) 



(4.4) 

(4.5) 
(4.6) 



where w(<p) — p/ p. The continuity equation can be interpreted as that of a fluid with a field dependent 
equation of state which in turns means a time dependent equation of state. 

We now present the equations for the linear perturbations in the Newtonian gauge for scales which 
are well inside the horizon following the notation of [65] . For a detailed derivation of the equations 
we refer to appendix (B). The system of equations governing the evolution of the linear perturbations 
is given as usual by the continuity equation, the Eulcr equation and the Poisson equation: 



^>«(w 



«j(r) S(k, t) + (1 + ew(r))kv{k, r) = 



v'(k,T) +H(T)v{k,T) + k<f>(k,T) = 

fc 2 *(fc, r) = 4TrGa 2 Q(k, r)p mc (r)(5(fc, r), 



(4.7) 

(4.8) 
(4.9) 



where 8 = Sp/ p and p mc is the minimally coupled background DM density. With this formalism the 
Q function takes the general form Q(k, r) = (1 + 0(e)), where the corrections come from the modified 
background DM density. 

From eq. (4.7) we can see that the continuity equation is modified in two ways: the last term can 
lead to a speed up or slow down of the growth of perturbations, depending on the sign of e; the second 
term on the left hand side is a new genuine effect of this model that closely resembles a dilution term. 

The Euler equation is modified in two parts: H{t) is modified as in the Friedman equation (4.1) 
and the gravitational potential is changed as from the Poisson equation. Notice also that despite the 
presence of an effective pressure, no Jeans length appears in the equation. This is a consequence of 
the time-like character of the pressure term, as noted above. Deviation from the ACDM model can 
be parametrized by two functions [66]: ( = (\f + <&)/4>, that characterizes the effects of anisotropic 



stresses, and Q related to deviation from the standard Poisson equation. The ACDM model has £ = 
and Q = 1. In our class of models there are no anisotropic stresses and hence £ = as well. However 
the Q function is in general different from unity and hence it is a measure of the departure from 
ACDM on which next generation surveys will cast strong constraints. As said, Q = 1 + ef(p,p) in 
general. Explicit forms for the function / will be given below. 

It is a standard procedure to take the derivative of eq.(4.7) and using eq.s (4.8) and (4.9) to 
obtain one single equation describing the evolution of the linear perturbations. In our case this gives: 



5"{k, r) - 3«'(r)e (j^f^ ~ w(r) j *(/,. , ) - :\H< 



(5p(k,T)\ , 



6(k,T) 



- 3H(T)e ( t P ^ ,T l - w(t)\ 6'(k, t) - ew'{r)kv{k, r) + (1 - ew(r))to/(k, r) = (4.10) 
\dp{k,T) J 

which reduces to the standard ACDM equation 6" = —kv for e — > 0. For completeness we give also 
the Euler equation in real space: 

v'(x, t) + H(t)v(x, t) = -V*(x, t). (4.11) 

Note again that here both H(t) and *&(x,t) are modified according to equations (4.1) and (4.9) 
respectively. 

5 Results 

In this section we present the results of the integration of equations (4.5)-(4.10) for different choices 
of the potential which directly translates in different time behaviors for the pressure. In the absence 
of a clear mechanism that can predict the form of the potential, we consider various examples and 
derive their possible cosmological implications. We will first explore the case in which the potential 
for the field (p is a constant and then we will consider a decaying exponential potential. These can 
be seen as the extrema of the general class of decaying potentials. In fact, all power law potentials 
would have an intermediate behavior, while a potential linear in the field would give an exponentially 
decaying time behavior. 

As a further work assumption we choose F = 1. This means that the coupling is only time 
dependent and not also scale dependent so that as soon as the critical temperature is reached the 
coupling is active everywhere. This is a crude simplification which enables us to illustrate, in a first 
concrete example, the impact that the dark matter non-minimal coupling illustrated in section (3) 
can have on observations. We leave more realistic scale dependent scenarios to future study. The set 
of functions and parameters used are reported in table ??. We also fix cosmological parameters as 
follows: D, A = 0.76, Q DM = 0.24 and H = 74 km/s/Mpc. 

Constant potential 

In the case of a constant potential the effects on the background density evolution are shown in figure 
(la) as compared to the evolution of standard cold dark matter. Remarkably the effect of the coupling 
is to mimic a cosmological constant as can be seen from the background solution: 

pW=pL(1 + ^) 3 + ^pL (5-1) 

P(z) = - e -^fp° dm (5-2) 

where e e // = £Pdm/P* w ith p\ m the present day density of dark matter, and p* a reference density 
characteristic of the scale under analysis. In this case the background dark matter fluid behaves as 
if it were composed by two fluids, one standard pressureless fluid plus a fluid with a cosmological 
constant equation of state p — —p, as pointed out also in [63] in a different context, so that there is no 
extra need to include a cosmological term in eq (4.1) The evolution for the density contrast is shown 



in figure (2a) as compared to that of the density contrast for the ACDM model. The suppression is 
enhanced for illustrative purposes but it is clear that the growth of linear perturbation is suppressed 
with respect to the one in standard ACDM model, in a way that resembles a ACDM model with a larger 
cosmological constant. This is also clear from figure (2b) where growth function / = — log 5/ log z for 
our model decreases faster compared to the ACDM function. 

The Q function, defined in the Poisson equation and responsible for deviation in the gravitational 
potential has the following form in the case of a constant potential: 

2(M=(i-^r) (5.3) 

which implies that the effective gravitational constant that generates the gravitational potential is 
reduced at those scales at which the interaction is active. 

To conclude, the effect of the gravitational self coupling with a constant dark matter potential is 
to add an extra contribution analogous to that given by a cosmological constant with the result that 
the growth of the density contrast is more suppressed than in the standard ACDM model. 

Exponential potential 

The exponential potential acts in a completely different way. The background DM density is shown 
in figure (lb) as compared to the evolution of standard pressureless dark matter. The effect of the 
non-minimal coupling is to slow down the dilution of the dark matter density as soon as the coupling 
is switched on. Notice that the sudden change in the evolution behavior is a mere consequence of the 
step function that switches on the coupling only for t > t c . More realistically, a smoother crossing 
is expected. The effect of the coupling fades away with time and the model tends asymptotically 
to ACDM . This is a welcome feature as it could possibly boost the number of high redshift clusters 
observed around z <~ 2 without affecting present day halos [67, 68], similarly to the scenario pictured 
in [69, 70] but here relying only on the dark matter non-minimal coupling to gravity. 
In this case the explicit form of the pressure term, as obtained from eq. (3.12) is: 



p(z) = ifarcsinh 
where 



Qa/^dm 
(l + z) 3 



4 / r, \ 2 



I /:!\' /3paV 



(5.4) 



^=2UJ » [H%) • <?*> 

In figure(3a) the evolution of the density contrast as a function of redshift is plotted . In this 
case the Q function has a complicated expression, not reported here, due to the non-trivial relation 
between time and redshift for a ACDM model. We just comment that also in this case the function is 
always less than one, thus reducing the gravitational potential. 

In figure (3b) we plot the growth function / as a function of redshift. As in the case of the 
constant potential, here again we notice that the effect of the coupling is to reduce the growth of linear 
perturbations. In this case, however, the sign of the coupling constant is opposite. The coupling e is 
a phenomenological parameter and different potentials subtend different theories thus the sign of the 
coupling is not a priori determined. A more mathematical explanation for this fact can be given in 
the limit of Einstein-de Sitter universe. In this case the equation to be solve for the density can be 
generally written as 

pi(z)-- z {p(z) + ep° M z a ). (5.6) 

This equation has the general solution 

p{z) = Po M z 3 + iCz a (5.7) 

where C = — e/(— 3 + a) with a ^ 3. For a constant potential a = and thus the particular 
solution to the differential equation is positive, hence the density is higher compared to ACDM . 

4 Notice that the case of a = 3 would simply rescale the coefficient of the homogeneous solution. 
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Figure 1: Background density plots for (a) the constant potential (left panel, red line) and (b) the 
exponential potential (right panel, red line), compared to standard pressureless dark matter (black 
dashed line). In both plots we fix F(X,ip) — 1 and z c = 5. The density is in g/cm 3 . The constant 
e e ff in the case of the constant potential is chosen in order to give a clear idea of the effects of the 
geometrical interaction term. In particular, in this case, e e // = 1. 
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Figure 2: (a) Evolution of the density contrast versus redshift in the case of constant potential (left 
panel, red line) and (b) the growth function f(z) = —dlogS/dlogz (right panel, red line), compared 
to the standard ACDM results (black dashed line), both with F(X, ip) = 1 and z c — 5. The constant 
e e ff in the case of the constant potential is chosen in order to give a clear idea of the effects of the 
geometrical interaction term. In particular, in this case, e e ff = 1- 

In the exponential case instead a — 6 and hence the situation is reversed. Interestingly, for an 
exponential potential, this suppression is limited in time: perturbations are maximally suppressed 
around the time of the switching but asymptotically the model reduces to ACDM . Again, This can 
be interesting when trying to get a higher number of halos in the past only, without affecting present 
abundances. 

6 Conclusions 

Several observational tests coming from future surveys will soon be able to further proof or disproof the 
ACDM model with respect to non-standard scenarios. In this paper we have addressed the question 
of whether the dark matter fluid can behave differently at galactic scales rather than at cosmological 
scales, due to the presence of a non-minimal interaction between dark matter and gravity. We have 
extended the analysis done in [57] and we have illustrated for the first time cosmological consequences 
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Figure 3: (a) Evolution of the density contrast versus redshift in the case of exponential potential, 
left panel red line and (b) the growth function f(z) = —dlogS/dlogz, right panel red line, compared 
with the standard ACDM results (black dashed line), both with F(X, <p) = 1 and z c = 5. 

of such a scenario, both at the background level and within linear perturbation theory. In particular we 
have shown that a non-minimally coupled dark matter fluid is able to produce two relevant effects: a 
pressure term for dark matter able to reduce the growth of structures at small scales, plus an effective 
interaction term between dark matter and baryons that can explain correlations between the two 
components of the cosmic fluid. 

In this scenario we have considered the situation in which dark matter, at suitably late times, 
undergoes some sort of phase transition (e.g. analogous to a Bose-Einstein condensation), consequently 
developing a coherence length of a size comparable to that of the local curvature radius, thus becoming 
non-minimally coupled. 

We have studied in details the dark matter pressure term, neglecting the roles of baryons in 
the present analysis. In particular we have analysed the system for two choices of the dark matter 
potential that generates the pressure term: a constant potential, resembling a cosmological constant 
contribution, and an exponential potential. These two choices are a good sample as all power potentials 
have intermediate behaviors. 

For a constant dark matter potential, the dark matter fluid behaves like the superposition of 
two fluids, one standard pressureless dust plus a fluid that behaves like a cosmological term, with a 
consequent suppression of the density contrast at small redshifts. 

In the case of an exponential dark matter potential the effects are mostly relevant near the 
time of activation of the coupling. The background density is enhanced and the linear growth is 
suppressed for a limited redshift interval. In fact the pressure decays with redshift so that standard 
ACDM evolution is recovered asymptotically. 

We also provided the Euler and Poisson equations (4.8), (4.9) in a form convenient for N-body 
simulations for any choice of the potential V together with explicit expressions and predictions for 
the parameters Q and £, that characterize the deviations from ACDM . In particular we found that 
C = 0, meaning that no anisotropic stresses are generated by our model, and Q < 1, after the coupling 
is switched on, thus reducing the gravitational potential. Constraints on these functions has been 
cast but they are rather weak and strongly model dependent [71-75] thus making necessary a direct 
confrontation between our model and observations in order to cast constraints on the deviations of 
our model from ACDM . 

For both choices of the potential we have obtained a suppression in the growth of linear per- 
turbations as in the figures (2) and (3). This is a good indication that this class of models may be 
a viable possibility to solve some of the ACDM paradigm problems, like the core-cusp, the missing 
satellites, high-z clusters. Of course, a non-linear analysis is required to evaluate these effects. 

We have derived the general perturbation equations valid for any F(ip,X) and V((p) (see eq.s 
(4.4, 4.5, 4.7, 4.8, 4.9)) , though we have limited our specific examples to F(cp, X) = I, for simplicity. 
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Baryons can be included, introducing a much welcome relation between dark matter properties and 
baryonic features as expected from observations [46]. 

We leave to follow-up studies a more realistic and quantitative description, including an environ- 
ment dependent coupling with a generic F((p, X); observations can be used to put constraints on the 
free parameters, e, p* and on the function F(ip,X). 
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A Derivation of the background equations 

In this appendix we give full details of the derivation of the equations reported in section (3). The 
starting action is 

S = ShbM + Sdm\$M + -Jd i x^XF(X,p) (x+ \v 2 (p)\ . (A.l) 

The variation with respect to the independent variables ip, A and the metric g^ respectively gives: 

ADv? + V M AVV - \V(<p)V v (<p) + — Uip (f(X, ip) Ux + \v(p)A + F V X (x + ^V 2 (p) 
V^W (f(X, p) Ux + \v{<pf\ + F(X, p) x X (x + X -V 2 {p) 

F(X, <p)XV(<p)V v {<p) + F(X, p) x X (x + l -V 2 {p>) 



0, (A.2) 



X--V(<p) 2 = 0, 

Guv = 8nGTu V + 8ttGT 



ltd 



where 



T, 



2e 



P = f t \ V " r-^rv 



F(X,p)i [ 2X+-V(p) 2 J + 



+ g^F(X 1 p)X[X + -V(p)' 



In order to simplify the notation we define the following quantity: 



1 

2" p* p 

so that the stress energy tensor can be rewritten as: 

T^ = (A + eCT c ) V^V^ + eC nmc g^ 



C nmc {X,p) = -F(X,p)T ^^ vLp = ffil^X ( X+ -V(p) 2 



(A.3) 
(A.4) 



(A.5) 



(A.6) 



(A.7) 



By a direct confrontation with the shape of the perfect fluid stress energy tensor we can identify the 
pressure of the field as: 

p = C nmc (X,p). (A.8) 
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Then, using the constraint equation (A. 3) we get: 

p = (A + eCT c ) V{ V f - eC nmc 7 (A.9) 

P = C nmc , (A.10) 

u' 1 = V(<p) _1 VV. (A.ll) 

The constraint (A. 3) can be rewritten as: 

<p=-V(<p), 's^V (A.12) 

while we have that: 

d = V^ = V(ip)- 1 Dcp + V v {ip) (A.13) 

with this set of definitions, equations (A. 2) and (A. 3) can be rewritten as: 

\ = V- 2 [V{ V )p^~(p + ep)d], (A.14) 

<p = -V(<p), (A.15) 

where we notice that the minus sign in equation (A.14) has appeared coming from the definition of 
derivative (A.12) and where: 

P , v = 2\vv ip + e [c n ^v 3 v v + cy^v 2 + £T c vv<p - £T C ] • ( A - 16 ) 

Here we have defined the various terms in the Lagrangian as follows: 

C— = ^,^, C!r = §*(*,„) ^ + l -F. x{ XM^, (A.17) 

2 p* 2 p* 2 p* 



p* 2 p* 2 p* 



CTx = 2 ± ^ + *F{X lV U x ^- + ±F(X )V ) x ^-, (A.18) 



^ = nx,^v M+ 1 -F vX ^L + 

P* & P* 

+ \F{X, v ) tX V tV (<p)Vy?- + \F(X, <P), V ^£, (A.19) 

<£ p* z, p* 

C nnic = ^X^Y^ftv^) + l F ( X ,<p) t Vyf. (A.20) 

£ p* £ p* 

Notice that the Klein-Gordon equation for the field tp is now an evolution equation for the Lagrangian 
multiplier A eq (A.14). 

B Linear cosmological perturbation theory 

In this appendix we derive the equations reported in section (4) following the notation of Kodama 
and Sasaki [65]. Given a metric g^ v , we perturb it around a background solution, 

9fiu = g^u + Sg^, (B.l) 

where Sg^ <C 1. The stress energy tensor of matter is perturbed in the same way around its perfect 
fluid form, 

f» = T» + 8T£ (B.2) 

To first order in the perturbed quantities the conservation equations V ' V T" = read: 

d^BT* + T'JTZ + ST^T? - T%8TS - 5T%T^ = 0. (B.3) 

In the Newtonian gauge the previous equations give for the time component 

6p'(r, k) + h (v(t, k)fc + 3*'(r, k)) + 3H(t)S p (t, k) + 2,M{t)5p{t, k) = 0, (B.4) 
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and 

h(T)v'(r, k) + (4W(r)/i(r) + ti{r)) v(t, k) - ekSp{T 7 k) - /i(r)fc*(r, k) + ^tt t (t, k)/c = 0, (B.5) 

for the spatial component. Perturbing the Einstein field equations (A. 4) one obtains a set of coupled 
equations that relates the perturbed metric quantities to the perturbed matter sources. The (0, 0) 
component describes how the perturbed matter density sources the gravitational potentials 

3-H(r) 2 *(r, k) - 3H(t)*'(t, k) - fc 2 *(r, k) = -4irGa 2 5p(T, k), (B.6) 

the (0, i) component relates gravitational potentials to the peculiar velocity of the matter, 

feW*(r, k) - fe*'(r, k) = 47rGa 2 /i(r)u(T, k), (B.7) 

while the trace of the (i, j) component relates gravitational potentials to the isotropic pressure, 

(2a" - H{t) 2 ) *(t, k) + H(t)*'(t, k) - \k 2 ^(r, k) - ^fc 2 *(r, k)+ 

- 2-H(t)*'(t, k) - *"(r, k) = 47rGa 2 e<5p(T, k). (B.8) 

As usual, the sum of the gravitational potentials is related to the anisotropic stresses: 

- k 2 (*(r, k) + *(r, k)) = 87rGa 2 ep(r)7r T (T, k) (B.9) 

where *f? is the Newtonian potential. Since the stress energy tensor is diagonal and the perturbed 
equations do not contain terms that mix different components, then ttt — and hence <& = — \&. 

Following [63] we can add another equation to the system, perturbing the constraint equation 
(A. 3). This is not independent as it turns out to be the equation for the velocity potential v = 
V((p)~ 1 k/ad(p. The constraint equation is 

{p' = g m u Q d (p, (B.10) 

where the tilde indicates the physical quantity containing also the perturbations. When perturbed in 
Newtonian gauge this gives: 

6tff = a\y v {(p)S(p + 9V{<p)]. (B.ll) 

Hence the system of equations to be solved is the following: 

6<ff = a(V v (<p)6tp + *V(<p)) (B.12) 

Sp' + h(kv- 3*') + ZH (Sp + eSp) = (B.13) 

H+-—}v-kV = ek— (B.14) 

2 p J p 

k 2 ^f + 3H(^ / + H^)=4TrGa 2 Sp (B.15) 

k (*' + %*) = 4irGa 2 hv (B.16) 

*" + H& + (ZH' + H 2 ) * = AirGa 2 eSp. (B.17) 

This set of equations is formally equivalent to the standard one apart from the first equation, derived 
from the perturbed constraint given by eq. (A. 3). However in the class of models under consideration 
extra relations exist that link together some of the variables. We have: 

v = V(ip)~ 1 k5<p, 

Sp = Pv S v . (B.18) 

Notice that the system (B.12)-(B.17) is redundant as we have 6 equations and 5 unknown variables. 
In what follows we rewrite Poisson and Euler equations (B.14), (B.15), in such a way that the modifi- 
cations due to the non-minimal coupling are manifest. 
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We are interested in the dynamics at scales well inside the horizon, A = H/k -C 1 and such that 
the characteristic evolution time for the potential & is of the order of T-L~ l . The system of equations 
governing the evolution of the perturbations in this limit reads: 

6(p' = a(V v (<p)5tp + ¥V(<p)) (B.19) 

dpi + hkv + iU (Sp + ep v S<p) = (B.20) 

H + e— ) v - k& = ek— (B.21) 

P J P 

fc 2 * = 4nGa 2 Sp (B.22) 

*' + W& = 4irGa 2 hv/k (B.23) 

*" + W + {pi! + "H 2 )* = AirGa 2 e8p. (B.24) 

Let's consider now Euler equation (B.21). Using (B.18) the two terms proportional to e cancel 
giving formally the same expression as in the standard case: 

v'(k, t) + U{T)v{k, t) - fc*(fc, t) = 0. (B.25) 

However here T-L is the modified Hubble expansion rate obtained by solving equation (4.1). Eq. (B.25) 
is therefore actually non standard both in the friction term and in the gravitational potential that 
feeds it, which, from eq. (B.22), can be rewritten as: 

fc 2 *(/c, t) = 4irGa 2 Q(k, T)p(r)6{k, r), (B.26) 

where the Q function measures deviation form the ACDM model which has Q = 1 and where 8 is the 
dimensionless density contrast 6 = Sp/p. Explicit shapes for this function are given in section (5). 

The continuity equation is usually rewritten in terms of the dimensionless density contrast S. 
This gives: 

S'(k, t) + 3eH(r) ( f^^l - w(t)] 5{k, r) + (1 + &vM) kv(k, r) = (B.27) 

\6p(k,T) J 

where w = p/p. The full analysis of these equations is carried out in section (3) and (4). 
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